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I.  INTRODUCTION 


This  report  describes  the  work  done  on  Contract  No.  N001 23-76-C-l 371 
for  the  period  1  February  1977  to  30  April  1977. 

The  overall  program  is  divided  into  three  areas.  These  are  1)  basic 
scattering  code  development,  2)  reflector  antenna  code  development  and 
3)  basic  theoretical  studies  to  support  the  first  two  areas.  The  following 
sections  describe  the  progress  made  during  the  second  quarter  in  each  of 
the  three  areas  mentioned  above. 


II.  PROGRAM  SCOPE 

The  scope  of  the  work  under  Contract  No.  N001 23-76-C-l 371  is  to 
develop  the  necessary  theory,  algorithms  and  computer  codes  for  simu¬ 
lating  antennas  at  UHF  and  above  in  a  complex  ship  environment.  A  mile¬ 
stone  chart  for  the  total  program,  which  extends  over  a  three  year  period, 
is  shown  in  Table  I.  A  more  detailed  breakdown  of  the  effort  planned  for 
the  first  year  is  shown  in  Table  II.  The  following  sections  describe  the 
progress  made  during  the  third  quarter  of  the  program  in  Table  II. 
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III.  BASIC  SCATTERING  CODE  DEVELOPMENT 


The  purpose  of  this  section  is  to  describe  the  present  status  of  the 
basic  scattering  code  development  for  the  analysis  of  antennas  in  a  complex 
shipboard  environment.  The  Geometrical  Theory  of  Diffraction  (GTD)  is 
being  used  to  develop  algorithms  to  solve  for  the  scattering  from  basic 
plate  and  cylinder  structures.  These  simple  components  can  then  be  com¬ 
bined  to  form  box-like  structures  with  nearby  finite  elliptic  cylinders 
that  can  represent  the  various  component  structures  of  a  ship.  The 
algorithms  are  being  implemented  into  a  user-oriented  computer  code. 

In  this  period  the  scattering  code  for  the  flat  plate  simulation  of 
convex  structures  has  been  completed  and  delivered  to  NOSC.  A  draft  copy 
of  a  user's  manual  also  was  delivered  and  the  scattering  code  has  been 
succesfully  implemented  on  the  NOSC  computer  system. 

The  early  delivery  of  the  computer  code  has  allowed  discussions  to 

take  place  that  will  aid  in  the  development  of  future  codes  to  better 
satisfy  the  specific  needs  of  NOSC  and  other  users.  For  example,  an  in¬ 
finite  ground  plane  option  has  been  added  to  the  code  now  under  develop¬ 

ment  based  on  the  needs  of  NOSC. 

The  methodology  used  to  input  the  data  into  the  scattering  code 
also  is  under  revision  at  present.  Default  input  data  are  being  placed 
into  the  program  for  specific  geometries.  Commands  are  being  supplied  so 
that  only  select  pieces  of  data  need  to  be  changed  for  any  particular 
computer  run. 

The  scattering  code  presently  under  development  will  incorporate  all 
of  the  above  improvements.  Also,  the  scope  and  generality  of  the  problems 
that  can  be  solved  will  be  increased.  The  code  will  allow  flat  plate 
simulation  of  convex  and  concave  structures,  finite  elliptic  cylinder 
structures,  and  a  psuedo  combination  of  plates  and  a  cylinder  structure. 
However,  the  coupling  between  the  plate  structures  and  cylinder  is  not 
included  at  this  time. 

The  scattering  code  for  the  flat  plate  simulation  of  concave  structures 
is  being  developed  to  allow  as  much  generality  as  possible  in  defining 
scattering  structures  by  flat  plates.  This  generality  also  leads  to  more 
complexity  in  the  number  of  coupling  terms  between  the  plates.  The 
dominent  scattered  field  components  up  to  the  second  order  of  interactions 
have  been  added.  This  would  include,  for  example,  double  reflected  fields 
between  the  plates,  etc.  If  a  particular  geometry  of  plates  is  desired 
that  requires  higher  order  scattering  terms  to  obtain  a  solution  useful  for 
engineering  purposes,  the  necessary  terms  can  be  easily  added  after  they 
have  been  put  into  the  GTD  format. 
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One  such  term  that  is  not  included  at  present,  but.  has  been  found  to 
be  important,  is  doubly-diffracted  fields.  A  theoretical  study  lias  been 
under  way  to  show  the  validity  of  the  fiTD  in  handling  double-diffraction. 
However,  if  the  doubly-di ffracted  fields  are  included  in  all  desired 
regions  of  space,  the  algorithm  would  be  very  time  consuming  and  hence 
costly  to  run.  This  can  be  overcome  by  writing  the  algorithm  in  such  a 
way  that  the  doubly-diffracted  fields  are  included  only  when  their  effects 
are  significant  in  the  antenna  pattern.  Techniaues  to  accomplish  this 
have  already  been  initiated,  but  to  fully  complete  and  test  them  will 
require  additional  tine  and  effort. 

The  scattering  code  for  the  finite  elliptic  cylinder  has  been  assembled 
and  is  presently  being  tested.  Improvements  were  made  in  the  algorithms 
this  past  period  to  improve  the  efficiency  and  extend  the  range  of 
applicability  of  the  code.  To  illustrate  the  wide  range  of  cylinder  sizes 
and  source  distances  that  can  be  calculated  with  this  scattering  code,  a 
series  of  radiation  patterns  for  the  geometry  shown  in  Figure  1  is  presented. 

The  first  set  of  patterns  is  for  a  X / 2  dipole  antenna,  mounted 
parallel  to  the  y-axis  at  a  fixed  distance  one  wavelength  from  the 
surface  of  a  perfectly  conducting  circular  cylinder,  the  antenna  pattern 
for  this  geometry  without  a  cylinder  present  is  shown  in  Figure  2.  With 
the  dipole  oriented  parallel  to  the  y-axis  the  energy  propagates  around  the 
cylinder  following  the  acoustically  hard  boundary  condition  where  the 
electric  field  is  perpendicular  to  the  cylinder  surface.  The  cylinder 
radius  is  increased  geometrically  by  a  factor  of  two  from  a  =  1/2X  to  64X 
and  the  results  are  shown  in  Figure  3.  All  the  results  are  shown  in  terms 
of  dB  plots  normalized  to  0  dB  with  10  dB  per  division. 

Note  that  there  is  a  slope  discontinuity  in  Figures  3f-h  at  approxi¬ 
mately  the  <P  =  90°  -  120'J  and  <p  =  240°  -  270°  regions.  This  is  due  to 
the  fact  that  the  dipole  field  can  produce  a  slope  field  associated  with 
the  transition  field  of  the  cylinder.  This  slope  transition  field  has 
not  been  added  at  this  time.  It  is  presently  under  development  for 
another  contract  and  can  be  added  when  completed  if  it  appears  necessary 
to  give  good  engineering  results. 

In  Figure  4  the  cylinder  has  a  fixed  radius  of  2X  and  the  source 
position  is  varied  1/2X,  2X ,  4X,  8X  off  the  cylinder  surface,  respectively. 

The  results  in  Figures  5  and  6  are  the  same  as  in  the  above  figures 
except  that  the  dipole  is  now  oriented  parallel  to  the  cylinder  axis. 

This  represents  the  acoustically  soft  boundary  condition  where  the  electric 
field  is  parallel  to  the  surface  of  the  cylinder.  Note  that  the  slope- 
transition  field  is  not  needed  in  this  case  because  the  source  field, 
which  is  isotropic  in  the  x-y  plane  for  this  case,  does  not  have  a 
slope  associated  with  it.  Further  tests  are  being  made  on  the  scattering 
code  and  they  will  be  presented  in  later  reports. 


The  psuedo  combination  of  plate  structures  and  a  finite  elliptic 
cylinder  is  also  being  tested  at  the  present  time.  This  part  of  the 
scattering  code  does  not  contain  the  interactions  between  the  plates 
and  cylinder;  however,  the  scattered  fields  from  the  plates  are 
shadowed  by  the  cylinder  and  vice  versa.  This  causes  discontinuities 
in  the  calculated  pattern,  but  this  can  be  very  infomati ve.  The  mag¬ 
nitude  of  the  discontinuity  gives  a  gauge  on  how  important  the  different 
coupling  fields  will  be  to  the  final  pattern.  A  start  has  been  made  on 
studying  some  of  these  coupling  terms  but  due  to  the  generality  the 
geometries  desired  and  the  complexity  of  finding  many  of  these  •  na 
great  care  will  have  to  be  taken  in  using  only  those  terms  nece  ry  for 
a  good  engineering  result  for  practical  structures.  Which  terms  -  use 
will  become  clearer  as  the  study  progresses  during  the  next  coni  t  yea 

A  user's  manual  will  be  completed  in  the  next  period  for  ti  ■  ates 
and  cylinder  scattering  code,  and  the  scattering  code  will  be  tesced  and 
prepared  for  delivery  to  NOSC  by  the  end  of  this  contract  year. 


— 4> 


d.  p  =  10X,  a  =  Z\ 


lounted  parallel  to  the 
iter  of  a  circular 
is  plotted  in  dB. 


IV.  REFLECTOR  ANTENNA  CODE  DEVELOPMENT 


The  purpose  of  the  present  effort  is  to  develop  a  user-oriented 
computer  program  package  by  which  the  far  field  pattern  of  a  typical  Navy 
reflector  antenna  can  be  calculated.  Feed  blockage  and  scattering  effects 
are  not  included  in  this  phase.  These  will  be  included  later  (see  Table  III). 

Development  of  the  applied  theory  to  be  used  in  the  computer  program 
package  for  the  general  reflector  antenna  has  continued.  As  reported  in 
the  Second  Quarterly  Report  (1  November  to  31  January  1977),  the  use  of 
overlapping  subapertures  with  triangular  distribution  provides  an  ac¬ 
curate  method  for  calculating  the  far  field  pattern  while  minimizing  the 
number  of  aperture  field  computations. 

The  computer  algorithms  for  the  aperture  integration  (AI)  part  of  the 
reflector  antenna  computer  code  for  the  far  field  pattern  are  being  programmed 
on  the  OSU  Datacraft  computer.  Programming  of  several  portions  of  the  com¬ 
puter  code  has  been  completed.  One  portion  involves  inputting  a  general  re¬ 
flector  rim  shape  and  then  subdividing  the  aperture  into  rectangular 
subapertures.  Another  portion  performs  the  2-D  numerical  integration 
required  to  calculate  the  far  field  pattern. 

The  theory  of  the  rectangular  subaperture  method  is  given  below. 

Following  that,  the  program  performance  specification  (PPS)  for  the 
reflector  antenna  computer  code  (FF  w/o  blockage)  is  given. 

A.  Rectangular  Subaperture  Method 

The  theoretical  approach  for  computing  the  far  field  pattern  of  the 
general  reflector  is  based  on  a  combination  of  the  Geometrical  Theory  of 
Diffraction  (GTD)  and  Aperture  Integration  (AI)  techniques.  AI  is  used 
to  compute  the  main  beam  and  near  sidelobes;  GTD  is  used  to  compute  the 
wide-angle  sidelobes  and  the  backlobes.  The  geometry  of  the  reflector 
rim  is  treated  as  piece-wise  linear;  that  is,  a  series  of  linear  seg¬ 
ments  is  used  to  approximate  a  smoothly  curved  rim  such  as  a  circular 
aperture. 

Although  GTD  is  used  to  calculate  most  of  the  far  field  pattern, 
the  main  beam  and  near  sidelobes  of  the  reflector  antenna  are  calculated 
by  the  aperture  integration  approach.  Since  the  required  integration  is 
usually  very  cumbersome  for  electrically  large  antennas,  an  approach 
based  on  using  overlapping  rectangular  subapertures  is  used.  Each  sub¬ 
aperture  can  be  electrically  large,  i.e.,  several  wavelengths  in  size. 

This  minimizes  the  number  of  subapertures  required.  The  theory  of  the 
rectangular  subaperture  method  for  aperture  integration  will  now  be  given. 

The  reflector  rim  geometry  is  specified  by  the  x-  and  y-coordinates 
of  each  junction  point  L|<  and  iin  on  the  lower  and  upper  portions  of  the  rim, 
respectively.  Thus  the  projection  of  the  reflector  rim  onto  the  aperture 
plane  is  specified  as  piece-wise  linear.  A  typical  reflector  rim  shape  is 
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shown  in  Figure  7a,  with  four  points  on  each  of  its  lower  and  upper  portions. 
A  rectangular  grid  size  (dx  and  dy)  is  then  chosen  so  that  the  aperture  can 
be  divided  into  rectangular  subapertures.  A  typical  subaperture  is  shown 
in  Figure  7a  where  the  aperture  field  distribution  is  specified  by  its 
value  at  the  center  of  the  subaperture  and  by  its  values  at  the  8  points 
on  the  edge  of  subaperture  (centers  for  the  overlapping  subapertures)  as 
shown  in  Figure  7a.  Thus  the  distribution  over  the  entire  aperture  plane 
is  specified  by  its  values  on  the  rectangular  grid  as  shown  in  Figure  7b. 

The  vertical  grid  lines  are  numbered  by  integer  values  of  I  =  1  to  Imax 
corresponding  to  the  x-coordinate  as  shown  in  Figure  7b.  The  horizontal 
grid  lines  corresponding  to  the  y-coordinate  are  numbered  by  the  integer 
J,  starting  with  J  =  1  for  the  lowest  point  on  the  rim.  The  center  of 
the  aperture  plane,  which  is  on  the  z  axis  of  the  reflector,  is  designated 
by  I  =  Ic  and  J  =  Jc  as  shown.  The  minimum  and  maximum  values  of  J  on 
each  vertical  grid  line  (I)  are  determined  by  the  rim  geometry. 

The  y-coordinates  of  the  intersection  of  the  lower  rim  and  the  vertical 
grid  line  (I)  is  designated  by  Yl  as  shown  in  Figure  8.  The  horizontal 
grid  line  closest  to  this  intersection  is  designated  by  J  =  Jlr  (I). 
Similarly,  the  horizontal  grid  closest  to  the  intersection  (at  Yu)  of  the 
same  vertical  grid  line  with  the  upper  rim  is  designated  by  J  =  Jur(I). 

In  the  computer  code  the  array  variables  Jlr(I)  and  J(jR ( 1 )  serve  as  the 
lower  and  upper  limits  of  the  y-integration  for  each  vertical  grid  line  (I). 

The  aperture  distribution  for  the  basic  subaperture  (Type  0)  is  a 
triangular  one  as  given  by 


(x,y) 


(1) 


for  |x-x0|  <  dx  and  |y-y0|  <  dy;  and  fSo  is  zero  otherwise.  The  resulting 
far  field  pattern  for  the  basic  subaperture,  i.e.,  the  element  pattern,  is 
given  by 


where 

4>x  =  k  d^  sin  0  cos4> 

and 

$y  =  k  d^  sin  0  sin<j>  . 
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The  phase  reference  for  Equation  (2)  is  taken  at  the  center  of  the  sub¬ 
aperture  (x0,y0);  the  phase  is  referred  to  the  origin  of  the  aperture 
plane  during  the  integration  process. 


UPPER  RIM 


Figure  8.  Lower  and  upper  horizontal  grid  lines  for  y-i ntegration. 
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TABLE  III 


Subaperture 
Type  N 


Interval 
for  x-x 

0 

Interval 
for  y-yQ 

Element  Pattern 

FSN  6 

(-dx ,+dx) 

(-dy ,+dy) 

dx  dy  FF  (4>x)  Fp  (<py) 

(-dx,+dx) 

(O.+dy) 

dx  dy  Fp  (<px)  FH  ($y) 

(-dx ,+dx) 

(-dy,0) 

dx  dy  Fp  (<f>x)  F^  (-<|>y) 

(O.+dx) 

(-dy ,+dy ) 

dx  dy  FH  ( +4> x )  Fp  (<py) 

(-dx ,0) 

(-dy ,+dy ) 

dx  dy  Fh  (-♦x)  Ff  ( 4>y ) 

In  order  to  treat  the  subapertures  along  the  reflector  rim,  four  other 
types  of  subapertures  are  introduced  as  shown  in  Figure  7b.  The  area  of 
each  additional  type  is  one-half  of  that  for  the  basic  subaperture  (type  0) 
The  field  distribution  for  each  type  is  the  same  as  that  given  by  Equation 
(1)  except  for  different  x  and  y  intervals  as  given  in  Table  III.  The 
corresponding  element  patterns  can  be  expressed  in  terms  of  the  one¬ 
dimensional  pattern  functions:  Ff  (4>s)  for  a  full  triangular  distribution, 
fp  (s),  and  Fp|  ( 4>s )  for  a  half  triangular  distribution,  fp|  (s),  as  shown 
in  Figure  9.  Thus 


Figure  9.  Triangular  subaperture  distributions  (full  and  half) 


The  resulting  element  patterns  F$n(6,<}>)  for  each  type  are  also  given 

in  Table  III . 

The  y  integration  is  divided  into  3  parts:  the  lower  part  usually 
consists  of  the  contribution  of  the  Type  1  subaperture  at  the  lower  rim 
and  the  upper  part  is  usually  one  subaperture  of  Type  2  at  the  upper 
rim.  Types  3  and  4  are  used  for  the  lower  and  upper  parts  when  the 
reflector  rim  has  a  steep  slope.  The  subapertures  in  the  middle  part 
are  the  basic  Type  0.  The  Type  3  subaperture  is  used  exclusively  along 
the  left  side  of  the  rim  (I  =  1);  and  the  Type  4  is  used  on  the  right 
side  (I  =  Imax),  as  shown  in  Figure  7b. 

Dividing  the  y  integration  into  3  parts  has  the  significant  advantage 
that  the  product  of  the  aperture  field  and  phase  exponential  for  each 
subaperture  can  be  simply  sunned  for  each  integation  part.  The  result  for 
each  part  of  the  y-integration  is  then  multiplied  by  the  element  pattern 
Fsn (®  »4> )  of  the  appropriate  Type  N.  Thus  the  contribution  from  each  part 
can  be  expressed  as 

JF 

Ys  =  FSN(e,*)  I  Ea  (lx  dx,  Jy  dy)  eJ  Jy  *y  (5) 

J=Jj 

where  lx  =  I - I c  and  Jy  =  J-Jc. 

The  exponential  in  Equation  (5)  provides  the  y  component  of  the  phase 
reference  for  each  subaperture.  The  lower  limit  J j  and  the  upper  limit 
Jp  are  given  in  Table  IV  for  each  integration  part,  where  N[_q  and  Njjq 


y-Integration 

Part 

lower 

^YSL^ 

Middle 

Upper 

^YSU^ 

TABLE  IV 


Subaperture 
Type  Number 


Limit 

Lower  (Jj)  Upper  (Jp) 


H  =  N, 


N  =  0 


N  NUP 
usually  2 


(usual ly  1 ) 

(Remaining) 

Jlr(D+n 

nug 

(usually  1) 

JUR( 

jlr<i>+nlg-1 


J,m  I  -N 


are  the  number  of  subapertures  in  the  lower  and  upper  parts,  respectively. 
The  values  of  N(_q  and  N(jq  are  usually  unity  except  for  very  steep  slopes 
of  the  reflector  rim,  where  the  subaperture  Types  N|_W  and  f4yp  may  be 
chosen  as  3  or  4  to  better  fit  the  aperture  rim  shape.  As  previously 
mentioned,  Type  3  is  used  exclusive’v  for  the  first  vertical  grid  line 
(I  =  1);  and  Type  4  is  used  for  the  final  grid  line  (I  =  Imax).  The 
y-i ntegration  sum  is  thus  given  by 


YSUM^  ” 


1 SL 


+  YSM  +  Y 


su 


(6) 


The  2-D  aperture  integration  is  completed  by  summing  on  I;  thus  the 
rectangular  subaperture  method  gives  the  far  field  pattern  as 


F(e,$) 


Imax 

ysum  ^ 


gj  (Ix)4>X 


(7) 


where  the  exponential  in  Equation  (7)  provides  the  x  component  of  the 
phase  reference.  Equation  (7)  represents  the  principally-polarized 
component  which  is  dominant  in  the  main  beam  region.  When  the  present 
far  field  code  is  delivered  to  NOSC  early  in  the  second  year's  program,  it 
will  include  the  cross-polarized  component  and  hence  the  total  field. 

6.  Program  Performance  Specification  ( PPS ) 

for  Reflector' Antenna  Code/ far  Field  w/o  Blockage 

The  objective  of  the  computer  code  is  to  provide  the  capability  for 
calculating  the  far  field  (FF)  patterns  for  typical  Navy  reflector  antennas. 
The  theoretical  approach  for  computing  the  far  field  pattern  of  the  general 
reflector  is  based  on  a  combination  of  the  Geometrical  Theory  of  Diffraction 
(GTD)  and  Aperture  Integration  (AI)  techniques.  AI  is  used  to  compute  the 
main  beam  and  near  sidelobes;  GTD  is  used  to  compute  the  wide-angle  side- 
lobes  and  the  backlobes. 

The  code  provides  the  following  capabilities: 

1.  Reflector  surfaces  may  be  either  paraboloidal  or  cylindrical- 
parabol ic. 

2.  A  general  reflector  rim  shape  may  be  used  (piece-wise  linear 
up  to  100  segments). 

3.  The  required  input  data  for  the  feed  pattern  is  minimized  by 
polynomial  and  piece-wise  linear  pattern  fitting. 

4.  Storage  and  computation  time  of  aperture  data  for  AI 
is  minimized  by  a  rectangular  subaperture  method. 
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5.  The  efficiency  of  the  far  field  pattern  confutation  is  max¬ 
imized  by  the  use  of  GTD  for  wide  pattern  angles  and  tne  use 
of  trie  rectangular  subaperture  method  for  near-axis  angles 
(main  beam  region). 

An  outline  of  the  code  is  given  in  Table  V.  Formats  for  input  and 
output  data  are  given  in  Tables  VI  and  VII,  respectively. 


TABLE  V 

OUTLINE  OF  COMPUTER  CODE 
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TABLE  VI 
FORMAT  OF  INPU1 


INPUT 

Choose  dimension  units  and  input 
all  linear  dimensions  in 
same  units 
i  .e. 

Focal  distance 
Rectangular  grid  size 


Dimension 


inches 

meters 

or  wavelengths 


Wavelength 


A  inches 
A  meters 
A  =  1 


Projected  rim  shape:  Coordinates  (Xn,Yn)  of  junction  between  each  piece- 

wise  linear  segment 


Parabola  Type:  Paraboloidal:  Type  1  (NREFL  =  1) 

Cylindrical:  Type  2  (NREFL  =  2) 


Input  all  angles  in  degrees: 


FF  Pattern  angle  data:  9,$ 


Feed  Pattern:  Pattern  level  in  dB: 

Angles  in  degrees:  ip,<p 
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TABLE  VII 
FORMAT  OF  OUTPUT 


All  input  data  for  information  and  checking. 

Complete  feed  pattern  as  calculated  from  input  feed  data. 
$,<!»  in  degrees;  f(<j>,4>)  in  dB 
Calculated  aperture  field  data 


Output  FF  Pattern 


<}>-angle  =  degrees 

Pattern  Angle 

e 

degrees 


FF  Pattern 

0  Component  <f> 

dB  Phase  dB 


Estimated  Machine  Storage  Required:  200  K  bytes 


Component 

Pha 


V.  THEORETICAL  STUDIES 


In  the  present  period  work  on  the  derivation  of  diffraction  coeffi¬ 
cients  for  a  perfectly-conductinn  circular  cylinder  of  radius  a  was 
completed.  This  diffraction  coefficient  is  particularly  useful  and 
numerically  efficient  as  ka  0;  furthermore,  it  may  be  used  when  the 
axis  of  the  cylinder  is  curved  as  in  the  case  of  a  bent  wire.  The 
dyadic  (or  matrix)  form  of  the  diffraction  coefficient  is  the  sane  as 
that  of  a  curved  wedge,  which  is  not  surprising  since  both  structures 
involve  lines  of  scattering  (the  axis  of  the  cylinder  and  the  edge  of  the 
wedge).  Although  diffraction  coefficients  were  obtained  only  for  cylinders 
of  circular  cross  section,  the  method  of  solution  can  be  readily  extended 
to  perfectly-conducting  cylinders  of  arbitrary  cross  section  and  to 
penetrable  cylinders,  provided  that  the  maximum  extent  of  the  cylinder 
cross  section  is  not  large  in  terns  of  a  wavelength. 

The  diffraction  coefficient  for  the  circular  cylinder  was  employed 
to  calculate  the  radiation  nattern  of  a  LAMPS  antenna  with  a  cylindrical 
mast  in  front  of  it.  A  separate  report  has  been  written  describing  this 
work  [1]. 

Work  on  the  GTD  analysis  of  the  scattering  by  two  staggered,  parallel 
half-planes  (plates)  referred  to  in  the  March  monthly  report  has  been 
extended  to  the  case  where  the  interacting  edges  form  a  thick  edge,  a 
geometry  more  likely  to  be  encountered  in  practice.  As  in  the  case  of 
the  staggered  half  planes,  the  second  edge  is  positioned  on  the  shadow 
boundary  of  the  first.  Some  difficulties  have  been  encountered  when  the 
second  edge  is  in  the  transition  region  of  the  first  edge  (but  not  on 
its  shadow  boundary),  and  the  field  point  lies  in  the  transition  region 
of  the  second  edge.  It  is  a  matter  of  reducing  a  double  integral  to 
a  form  suitable  for  numerical  results.  Also  in  the  present  period  we 
were  concerned  about  the  complex  angles  of  incidence  which  occur  in  our 
integral  representations  of  the  diffracted  field.  In  particular,  could 
we  justify  the  use  of  diffraction  coefficients  with  complex  angles  of 
incidence?  We  have  examined  this  point  carefully,  and  modified  our 
approach  so  that  this  difficulty  appears  to  be  overcome. 

In  the  period  ahead  we  plan  to  return  to  the  task  of  completing 
the  report  on  slope  diffraction. 
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